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Finite temperature spectral function of a hole in a quantum antiferromagnet and role 
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We study thermal broadening of the hole spectral function of the two-dimensional t — J model 
(and its extensions) within the non-crossing approximation with and without the contribution of 
optical phonons. We find that phonons at finite temperature broaden the lowest energy quasiparticle 
peak, however, the string excitations survive even for relatively strong electron-phonon coupling. 
Experimental angle resolved photo-emission spectroscopy (ARPES) results compare well with our 
calculations at finite temperature when we use strong electron-phonon coupling without any adhoc 
broadening. In addition, we have studied the role of vertex corrections and we find that their 
contribution allows us achieve the same overall agreement with the ARPES experimental results 
but using smaller values for the electron-phonon coupling. 

PACS numbers: 71.10.-w,71.10.Fd,71.27.+a,74.72.-h,79.60.-i 



I. INTRODUCTION 

The cuprous oxide superconductors show a broad peak 
near the Fermi energy followed by a "waterfall" -like fea- 
ture at higher energies in rather recent high resolution 
angle-resolved-photoemission-spectroscopy (ARPES) 
measurement^i^i^. Calculations based on the t — J 
model give a well-defined quasiparticle-like low energy 
peak and higher energy "string-like" excitations^^^. The 
results obtained from the t — J and the t — t' — t" — J 
models, using an artificial broadening of the lowest 
energy peak and of the other peaks corresponding to 
the string-excitations, agree reasonably well with the 
experimental spectra^. Furthermore, there are similar 
studies using the Hubbard model and its extensions^i^ii^ 
also indicating that the above features seen in the 
ARPES studies could be due to higher energy hole 
excitations arising naturally in these strongly correlated 
electronic models. 

In this paper we consider the role of finite temperature 
and of the coupling of the hole to optical phonons, as re- 
cent experiments have provided increasing evidence that 
electron-phonon coupling is strong in cuprates^^^^^^^^^. 
Our goal is to examine (a) whether or not the string 
excitations, claimed in Ref. i6| to be the cause of the 
"waterfall" -like features seen in the ARPES studies, sur- 
vive the presence of such strongly coupled phonons and 
(b) whether or not a natural broadening mechanism 
due to (i) finite-temperature and/or (ii) the coupling to 
phonons can give a reasonable explanation of the ob- 
served features of the ARPES spectra. 

Calculations based on the t — J model at finite tem- 
perature have been done using the Lanczos method^^, 
quantum Monte Carlo (QMC)^^ and recently using the 
so-called hybrid dynamical momentum average (HDMA) 
method^^. While the results obtained from these meth- 
ods are quite useful, the conclusions drawn from any one 
of them should be taken with some caution; for exam- 
ple, the Lanczos method can be applied to very small 



size lattices, and the so-called maximum-entropy tech- 
nique which is utilized by QMC disregards the high en- 
ergy peaks due to string excitations and other important 
details of the spectral function. Though the recently used 
HDMA method takes into account the electron-phonon 
vertex diagrams and produces the broad lowest energy 
peak quite well, the values of the coupling constant 7 
considered is not strong enough to justify the applica- 
tion of the momentum average method^^. 

In this paper we extend the method introduced in 
Ref. y and developed in Ref. 5ll6lll3l, at finite temper- 
ature and we also include the role of optical phonons. 
In Refs. IjIH the boson degrees of freedom were treated 
within the so-called spin-wave approximation and their 
coupling to electron and hole degrees of freedom is lin- 
earized with respect to boson creation operators. Fur- 
thermore, the self-consistent Dyson's equation for the 
single-hole spectral function was solved within the so- 
called non-crossing approximation (NCA) where only 
topologically "planar" diagrams are retained. In the 
present paper we work within the same linearized Hamil- 
tonian and we include the linear coupling to optical 
phonons as captured by the Holstein electron-phonon in- 
teraction. The calculations are carried out at finite tem- 
perature by solving the Dyson's equation within the NCA 
for both diagrams which include propagation of spin- 
wave excitations and diagrams which include phonon 
propagation. We find that together the phonons with the 
inclusion of the thermal broadening at room temperature 
give rise to a broadened spectral function which exhibits 
similar characteristics to those found in the ARPES stud- 
ies. More precisely, the conclusions of Ref. @ are valid 
without the need to artificially broaden the spectral 
function. Furthermore, as it is well-known the leading 
vertex correction due to coupling to spin-waves is zero 
and other higher order vertex corrections give negligi- 
ble contribution^^. In the present calculation we include 
the leading (two-loop) vertex corrections due to the hole- 
phonon coupling and due to the coupling of the hole to 



spin-wave excitations and we conclude that their contri- 
bution ahows us achieve the same quahtative agreement 
with the experimentahy determined hole spectral func- 
tion using smaller values of the electron-phonon coupling 
constant. 

In the following section (Sec. [Ill) ^^ describe the for- 
malism and the approach. In Sec. IIIII we present our 
results for the spectral function obtained by a numerical 
solution of the Dyson's equation. In Sec.[lVlwe compare 
our results with the experimentally determined spectral 
function. In Sec. [Vl we include the contribution of the 
vertex corrections and in Sec. [Vll we present the main 
conclusions drawn from the present study. 



II. FORMULATION 

The motion of a single hole in a spin-^ Heisen- 
berg antiferromagnet^ in a two-dimensional(2D) square 
lattice has been extensively studied using the two- 
dimensional (2D) t — J model^: 



J. 



Y^[js^^s^, + ^{sts- + s-sj)]. 



(1) 



<«,i> 



The first term is the usual hole-hopping term which op- 
erates in a space of singly occupied sites and the sec- 
ond and third form the usual Hamiltonian of the Heisen- 
berg antiferromagnet where we have allowed for a pos- 
sible anisotropy of the coupling between the z and the 
perpendicular spin components. While this Hamiltonian 
has been thoroughly studied during the last almost two 
decades using many techniques, only a handful of meth- 
ods are shown to yield accurate results in certain lim- 
its. For the case of a single hole spectral function, one 
of such rather successful techniques is the so-called self- 
consistent Born approximation (SCBA)^. 

A simple way to introduce the coupling of the hole 
motion to a single optical phonon branch is by adding 
to the t — J model an electron-phonon coupling term by 
means of the following Holstein term: 



Hel-ph = 1^0 ^ h\hk + -7= ^ c\ck-qhq + h.c. 

k ^ k,q 



(2) 



where b^ is the optical-phonon creation operator, Qq is 
a characteristic optical phonon frequency and 7 is the 
electron-phonon coupling constant. 

Within the linear spin-wave approximation, using the 
Bogoliubov transformation to diagonalize the Heisenberg 
term and and by linearizing the hopping term with re- 
spect to the spin-deviation operators, one finds^*^ the 
following expression for the Hamiltonian given by Eq. [T} 

H = ^o + J^(/,l/k + 4^k)+^t^k(4«k + /?,^A) 

k k 



+ Yhlfk-q[g{k,q)aq^g{k-q,-q)plq\ 

k,q 

+ flhk-q[g{k-q,-q)alq^g{k,q)/3q]^H.c. (3) 

The function g{k^q) which plays the role of the hole-spin 
wave coupling constant is defined in Refs. [3|5|. 

In order to calculate the effects of finite temperature 
we will use the Matsubara technique followed by ana- 
lytic continuation to the real frequencies to obtain the re- 
duced Green's function^^^^^. The self-consistent solution 
to Dyson's equation for the self energy of the Hamilto- 
nian given by Eq.[3]is obtained by iterating the following 
equation with respect to n: 

S(^+i)(k,cj) = ^^^(k,q) ArqGn(k-q,cj + cjq) + 



(1 + Arq)Gn(k-q,Cc;-cjq) + 

^ de 

—nF{e)Do{q, e - uj)ImGr(k - q, el\L) 



Gn(k,cj) 



1 



CJ - ^k - Sn(k,Cj) 



(5) 



where A^q = ^^(cc^q), with nsiu;) = ^^J_-^ , and nF(^q) = 
f^^l,-, ' Also, (,k = ^k — M, where e^ is the zeroth order 
hole energy, which according to Eq. [3] is equal to J, and 
/i is the chemical potential. Here oOq is the spin-wave fre- 
quency, 1^0 is the spin- wave propagator given as follows: 



Do{^,Lo)- ^^' 



CO'^ — COq 



(6) 



and Gr{uj^ k) is the retarded Green's function which is ob- 
tained from the Matsubara Green's function by analytic 
continuation: 



Gr(k, cj) = lim Giuj - 
5^0 



i5^ k). 



(7) 



In the case where we include the coupling to the opti- 
cal phonons via Eq. [2] we need to add to the above ex- 
pression for the self-energy three more terms which are 
the same as the above and they are obtained from the 
above expressions by replacing the hole- spin- wave cou- 
pling constant g by 7/a/]V, and the spin-wave frequency 
by the phonon frequency VLq. 

While the first two terms of the above equation are of 
order unity, it can be shown that the last term is of order 
of 1/A^ because we only consider a single hole. This can 
be seen by considering the following identity 



/CO 
nF{oj)A{p,uj). 
-CO 



(8) 



The difference between this equation and the third term 
of Eq.[4]is the presence of boson propagator, Dq which is 
not an extensive quantity. Hence, the order of magnitude 
of the last term in Eq. [4]is ~ 0{Nh/N) and vanishes in 
the thermodynamic limit. 
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FIG. 1: The calculated spectral function for k = (^, ^) fc 
pt = 10 for a 4 X 4(top), 8 x 8(middle) and 16 x 16 (botton 
lattice with (dashed line) and without (solid line) the inclusio 
of the third term in the self-energy expression given by Eq. (|4] 



The vanishing of the third term in Eq.[4]is also demor 
strated numerically in Fig. [1] which shows the spectrj 
function for (f , f ) with and without the third term i 
Eq. Hfor pt = 10 and for 4 x 4, 8 x 8 and 16 x 16 lattices 
respectively. The solid lines are the spectral functions 
without the third term while the dashed lines are ob- 
tained by including it by means of a single iteration of 
Eqs. 14151 Notice that for large enough size lattice the 
contribution of this term becomes negligible. In the rest 
of our calculations presented in this paper this term will 
be neglected. 

The self-consistent Dyson's equation in conjunction 
with the so-called non-crossing approximation(NCA) 
(crossing diagrams have a small contribution as explained 
in Refs. .511131) is solved by means of an iterative proce- 



dure to obtain the dressed hole propagator and the hole 
spectral function. For numerical calculations a small con- 
verging parameter 77 is needed in the zeroth order Green's 
function as follows 



G^''\k,co) = 



u; -^k^ir] 



(9) 



Starting from the above zeroth order approximation for 
the single hole Green's function, the Dyson's equation 
is iterated until convergence is achieved. Because the 
lowest-energy quasiparticle peak corresponds to a well- 
defined excitation, its width and height, as smaller and 
smaller values of r] are used, scale proportionally to r] and 
l/f] respectively. In addition, in order to avoid finite-size 
effects a smaller value of r] requires a larger size lattice. In 
Ref. 5 it was demonstrated that the single hole spectral 
function has negligible finite-size effects for lattices larger 
than 16 x 16 when a value for rj = O.lt was used. However, 
when we take smaller values of 77, we need bigger size 
lattices to eliminate the finite-size effects. For example, 
at T = and without phonons, we have found that in 
order to reach the thermodynamic limit for r]/t = 0.1, 
a 16 X 16 size-lattice is large enough, while if we take 
r]/t = 0.05 or r]/t = 0.01 lattices of sizes 24 x 24 and 
32 X 32 respectively are required. 
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r|/t=0.005 
r| /t = 0.01 




FIG. 2: The single- hole spectral function at T = with 7 = 
for (|, I) calculated on a 64x64 lattice with 77 = 0.005t (solid 
line) and O.Olt (dashed line). 

Since the quasiparticle peak is a well-defined 
excitation^ at T = 0, its height becomes greater as we 
decrease rj. On the other hand in the higher energy part 
of the spectral function, because it forms a continuum 
of states, the various Lorentzian contributions overlap 
strongly because the energy spacing Ae of neighboring 
energy levels becomes exponentially small with lattice 
size, i.e., Ae ~ g-^^ (where N = L x L is the number 
of lattice sites). Fig|2] shows the (f , f ) spectral func- 
tion with T] = O.Olt and 0.005t for a 64 x 64 size lattice. 



The difference is mainly in the height of the lowest en- 
ergy peak which doubles by decreasing the value of rj/t 
by a factor of 2 and the other parts of both of the spec- 
tral functions are very close. In this paper we have used 
J = 0.3t and ry = 0.05t. 



III. NUMERICAL RESULTS 

A. Finite temperature no phonons 

First of all we study the effect of temperature alone, 
i.e., without any phonons in the system. Fig. [3] 
shows the spectral function for k = (f,f) for T/t = 
0.001,0.01,0.05,0.1 and 0.15, calculated on a 16 x 16 
lattice. As the thermal broadening is most prominent 
near the lowest energy well-defined peak, notice that the 
multi-peak structure of the spectral function just above 
the lowest peak becomes more and more broadened as the 
temperature is raised. The effect of finite temperature is 
also to move the low energy peaks towards lower energies 
(a shift of about 0.05t occurs for T = O.lt). In Fig. H 
we present the spectral function for k = (0,0) using the 
same values of /3. Notice that the peaks which correspond 
to string excitations are robust even for temperature as 
high as T = 0.15t for both cases of the spectral functions. 
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FIG. 4: Spectral function for k = (0, 0) without phonons for 
a 16 X 16 lattice using pt = 1000, 100, 20, 10 and 6.67. 



and make the intensity graph similar to those observed 
in ARPES. 



B. Optical Phonons at T = 




FIG. 3: Hole spectral function for (^, ^) without phonons for 
a 16 X 16 lattice and for pt = 1000, 100, 20, 10, and 6.67. 

An intensity plot for /3t = 10 is presented in Fig. [5] 
along with the ARPES intensity^. Notice that there is a 
significant gap or pseudo-gap between the lowest energy 
peak and the first string excitation and also between the 
first string excitation and the peak which evolves to be- 
come an intense peak near (0, 0). In the following we will 
discuss that the presence of optical phonons which cou- 
ple strongly to the hole excitations can remove these gaps 



As it has been demonstrated the NCA is a good ap- 
proximation for the case of the single-hole in the pure 
t — J model, we need to make sure that it also works for 
the case of the perturbative expansion involving terms 
in which the boson loops are due to the phonon propa- 
gator. Mishchenko et al.^^ have shown that by increas- 
ing the value of the electron-phonon coupling strength 7 
and at zero temperature, a cross-over between the lowest 
state and the next string state takes place at 7 ~ OAt 
and from there on the lowest state which is like a narrow 
quasiparticle peak always stays dispersionless. Accord- 
ing to their calculation the next high energy state shows 
broadening just as it also appears in the experimental 
ARPES plot^. Surprisingly the results of Mishchenko 
et al.^^ with no phonons do not completely agree with 
earlier numerical results^. Here, we will study the ef- 
fect of non-crossing diagrams at finite temperature using 
this t — J-Holstein model with 7 both below and above 
the cross-over point ^^ and then compare our results with 
experimentally obtained intensity plots. 

Experimental values of characteristic phonon energy 
scales vary from 30 — SOmeV^'^^. It is rather well known 
that the value oft is approximately OAeV for the cuprate 
materials and we have used fto ^ OAt and 0.2t both 
producing essentially the same spectral functions^^. 

In Fig. [6] the spectral function for k = (f , f ) and for 
Qo = OAt with 7 = 0.0,0.1,0.2,0.4,0.5 and 1.0 calcu- 
lated on a 16 X 16 size lattice is presented. These results 
demonstrate that the string excitations are quite robust 
in the presence of optical phonons. In addition, in Fig. [3 
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FIG. 5: Top: Calculated intensity plot on a 48 x 48 lattice 
at /3t = 10 and no phonons. Bottom: ARPES intensity along 
the (0,0) to (7r/2,7r/2) direction. 



the k = (0,0) spectral function is shown, for Qq = O.lt 
with 7 = 0.0,0.1,0.2,0.4,0.5 and 1.0 as calculated on a 
16 X 16 size lattice. The same conclusion about the ro- 
bust nature of the string excitations can also be drawn 
from this graph. 

The spectral function calculated for J = 0.3t has some 
saw-teeth-like features at energy just above the lowest 
peak. Particularly we can see that the 1st such peak clos- 
est to the lowest energy quasiparticle peak (~ 0.2t energy 
apart) gains weight as the value of 7 is increased more 
and more, a feature which is also observed in DMC^^^^^ 
and in the HDMA^^ calculations at zero temperature. 
By increasing the value of 7 the low energy peaks move 
towards lower energies. For values of 7 of about O.bt this 
energy shift becomes approximately O.lt. When 7 = l.Ot, 
however, this energy shift becomes large (~ 0.4t). How- 
ever, NCA is not expected to be a good approximation 
to describe the spectral function for 7 = l.Ot. Also NCA 
does not show the cross-over phenomena as referenced by 
Mishchenko et al.^^. 

Phonons cannot broaden the lowest peak at zero 
temperature (due to energy conservation requirement) 
though the higher energy peaks corresponding to string 
states are broadened more and more with 7. In the next 
subsection (and Fig. [8|) we discuss that thermal broad- 




FIG. 6: Zero-temperature spectral function at k = (f , f ) on 
a 16 X 16 size lattice with optical phonons and Qo — O.lt and 
for 7 ranging from to 1.0. 




FIG. 7: Calculated spectral function for k = (O5 0) and 
at zero temperature for a 16 x 16 size-lattice for 7 = 



ening plays an important role in filling up the gap be- 
tween the lowest peak and the next phonon-generated 
small peak making an overall broad lowest energy peak. 



C. Finite temperature and optical phonons 



Fig. [8] shows the vicinity of the lowest energy peak 
of the calculated spectral function for k = (f,f) for 
P = 1000,100,20,10 and 6.67, and with ^0 = O.lt and 
7 = O.bt (Fig. Ela) and 7 = l.Ot (Fig. Elb) on a 16 x 16 
lattice. Notice that for both cases of electron-phonon 



coupling, as the temperature rises the lowest peak and 
the next phonon-induced peak smears due to thermal 
broadening and this gives rise to a single broad peak as 
seen in Fig. [HI This is the mechanism by means of which 
the lowest energy peak acquires a width of the same size 
as that found experimentally. A width of similar mag- 
nitude was used in Ref. [y to obtain agreement with the 
ARPES intensity. 
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FIG. 8: The spectral function for k = (f 5 f ) calculated on 
a 16 X 16 size lattice for pt = 1000, 100, 20, 10 and 6.67 at 
7 = 0.5t(top) and 7 = l.Ot (bottom). Only the vicinity of the 
lowest energy peak is shown. 

The effects of the electron-phonon interaction are pre- 
sented in Fig. [HI in a much wider frequency range, for k = 
(f,f) and for pt = 10 and for 7/t = 0,0.1,0.2,0.4,0.5 
and 1.0, where Qq = O.lt was used in this calculation. 
In Fig. [To] the spectral function for /c = and electron- 
phonon coupling 7 = 0.5t and for various values of tem- 
perature is shown. 

In Fig. [TTlf top) we present the calculated dispersion of 
the lowest energy quasiparticle peak for the t — J-Holstein 
model for various values of the electron-phonon coupling 
in the range j/t = — 1. In Fig. [TTlf bottom) we present 
the same calculation carried out for the t — t' — t" — J 
model for the parameter values believed to be needed in 
order to reproduce the ARPES data^^^^. Notice that the 
only significant effect on the t — t' — t" — J model, is to 
shift the overall energy by a constant and does not alter 
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FIG. 9: The spectral function for k -= (f , f ) for /3 = lOt on 
a 16 X 16 size lattice for 7 = 0.0, 0.1, 0.2, 0.4, 0.5, and 1.0. 
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FIG. 10: The spectral function for k = (0, 0) for 7 = 0.5t cal- 
culated on a 16 X 16 size lattice and for pt = 1000, 100, 20, 10 
and 6.67. 



the features of the dispersion. The effect of phonons on 
the hole dispersion for the case of the pure t — J model 
is more significant (Fig. [TT]). 



IV. COMPARISON WITH ARPES 

The observed ARPES spectral function reveals that 
the lowest energy peak for k = (f,f) has a width 
~ OAeV and it is the most intense feature together with 
the one near k = (0, 0). As we move from (f , f ) towards 
(0,0), the lowest energy peak moves gradually towards 
higher energies and a second peak grows near k = (0, 0). 
At around k = (^, ^) the intensity of the lowest energy 




FIG. 11: The hole band for f3t = 10 and for various values of 
7 for (top) the t — J model and (bottom) the t — t' — t" — J 
model with t' = -0.33t, t" = 0.22t. The values 7/t = 
0.0,0.1,0.2,0.3,0.4,0.5,0.65,0.8 and 1.0 were used in both 
graphs 



peak decreases appreciably and also smears over a region 
of higher energy. It is possible to explain these observa- 
tions using the results of the calculation based on the t—J 
model and the NCA as reported in Ref. 6 where the con- 
tribution of the string excitations gives rise to rather well- 
defined peaks in the spectral function at higher energy, 
provided that these string excitation peaks broaden sig- 
nificantly at around (f , f ) to give rise to some rather flat- 
intensity region. Furthermore, near k = (0,0), the peak 
which corresponds to a higher energy string excitation 
suddenly picks up intensity while the broadening process 
of the other string excitation peaks still prevails. This 
combined process of spectral- weight transfer and broad- 
ening of the peaks gives rise to the observed energy kinks 
in the ARPES intensity. However, these have been re- 
ported to be due to the electron-phonon interactions^^^^^ 
and the two energy scales separating the intermediate 
smeared intensity region from the two peaks on the two 
sides (one at (^, ^) and the other at (0,0)) has been iden- 
tified as the threshold of disintegration of the low-energy 
quasi-particles into a spinon and a holon branch^. 

The role of temperature is to broaden the high en- 
ergy string excitation peaks though the effect is not 
very pronounced^^ without simultaneously introducing 
the electron-phonon coupling. Fig. [T2l shows the spectral 




FIG. 12: The top-left, top-right and bottom-left spectral func- 
tions are along (0,0)^ (f 5 f ) ^^^ they are calculated for 
electron-phonon coupling 7/t = 0.2,0.5 and 1.0 respectively 
on 48 X 48 lattice. The bottom-right spectral function is the 
one obtained experimentally in ARPES. 



function for all values of k along the (0,0) -^ (f ? f ) di- 
rection on a 48 X 48 lattice for jS = 10, Qq = O.lt and for 
electron-phonon coupling j/t = 0.2, 0.5, and 1.0 and it is 
compared with experimentally obtained ARPES spectral 
function^. For k around (^, ^), where the height of the 
lowest peak and the higher energy string excitations are 
almost equal, the combined effect of the electron-phonon 
coupling and the finite temperature produces flat regions 
with much less intensity (than that of the lowest peak) 
in the spectral functions. Notice that using 7 = 0.2 and 
0.5t, which are below and near the possible cross-over 
point ^^, the high energy peak structure due to string ex- 
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FIG. 13: Intensity plot of a 48 x 48 lattice at /3 = 10 for 
7 = 0.2 (top) and 0.5 (bottom) 
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FIG. 14: The theoretical intensity plot for a 48 x 48 lattice 
at /3 = 10 and for 7 = 1.0 (top) is compared with the experi- 
mentally obtained ARPES spectra (bottom). 
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FIG. 15: Intensity plot calculated on a 48 x 48 lattice for 
7 = 0.5t and fit = 6.67 
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FIG. 16: Intensity plot calculated on a 48 x 48 lattice at 
7 = 0.5t and for pt = 10 with t' = -0.33t, t'' = 0.22t 



citations becomes visible at around k = (^, ^) and it is 
broadened. This mechanism creates a flat, low intensity 
region which becomes more pronounced for 7 = l.Ot. We 
note that in order for our NCA based calculation without 
vertex corrections to produce spectral functions and in- 
tensity plot (Fig. [Hj) similar to that obtained by ARPES 
we need to increase the value of 7/t to 1.0 (for pt = 10, 
which corresponds approximately to room-temperature). 
It also can be observed that as we move from (f , f ) to 
(0, 0) the peak structure which corresponds to string exci- 
tations becomes compressed more and more as a function 
of energy and the peaks are closest to each other at (0, 0). 
The same features can also be seen in Figs. [TSUH which 
present the intensity plot on a 48 x 48 lattice (for P = 10, 
Qq = O.lt) and for the same three values of 7/t = 0.2 (top 
part of Fig. [I3l), 0.5 (bottom part of Fig. [13]), and 1.0 
(top part of Fig. [14]) and it is compared with experimen- 
tally obtained ARPES spectral functions (bottom part 
of Fig. [Ej). Notice that the intensity plot becomes com- 
parable to the experimentally obtained ARPES spectral 
function when we use a large value of the electron-phonon 
7. As discussed in the following section, the presence 
of vertex corrections due to electron-phonon-coupling al- 
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lows us to use a smaller value of 7 to achieve the same 
qualitative agreement with the ARPES intensity. 

Since the value of t is not accurately known, we do not 
have a precise knowledge of the value of room tempera- 
ture in units of t. Hence, we also tried higher tempera- 
ture, T = 0.15t, where more broadening is obtained. As 
can be noticed from Fig. [151 where j/t = 0.5 was used, 
temperature could be an additional factor which helps 
us achieve better agreement with the observed intensity 
plot without having to increase the value of 7 and enter a 
domain where the validity of NCA becomes questionable. 

As we are mainly interested for intensity plots along 
the (0, 0) to (7r/2, 7r/2) cut in the k-space, the inclusion of 
t' and t'^ also does not change the spectral functions very 
much, as can be seen in Fig. [HI where the t — t' — t" — J 
model was used, taking t' = — 0.33t and t" = 0.22t and 
in addition /3 = 10 and 7 = 0.5t. 



culations indicating the discrepancy between the spectral 
function with and without vertex correction in the strong 
phonon coupling regime at zero temperaturei^*^ . In this 
section, we present the results of our study of the role 
of such vertex corrections, namely we improve the NCA 
by including the leading-order vertex corrections due to 
the electron-phonon coupling. The contributions to the 



gi:ci2 




V. VERTEX CORRECTIONS 

The electron-phonon vertex corrections are expected 
to be important for the hole-spectra unlike the electron- 
spin-wave vertex corrections whose contribution have 
been found to be small^-. There are some recent cal- 



FIG. 17: The leading vertex correction to the hole Green's 
function due to its coupling to phonons. 



self-energy are shown in the diagrams of Figs. llTlTSl and 
following the procedure outlined in Ref. [20 we obtain the 
following expressions: 



J 



S(")(cj,k) 



cr=±l qi,q2 



qi) 



+ (l + <^qJG(c.-ac.£_^^ 



N. 



(a) 



Xp) 



,(«) 



(a) 



G{u - a<4^ + <^_^^, q2)G(c. + <^_,,,, k - qi + q^) 



,(«) 



^,q2)G(cJ-cj, 
I 



(a) 
qi-q2- 



k-qi+q2)l/("Hk,qi,q2)A^!^q,. (10) 



The index a = 1,2 is used in order to distinguish the 
two different self-energy diagrams depicted in Figs. 1171181 
The two different cases of a are obtained as follows: 

• a=l. For the diagram depicted in Fig. [17] which 



involves only phonon loops 



,M 



,b) 



where 

Qo is the phonon frequency which we take it 
to be a constant characteristic optical phonon fre- 
quency Qq. In this case /^"^(k, qi,q2) = j^/N'^. 



a=2. For each of the diagrams depicted 

in Fig. [18] which involve one phonon and 
one spin-wave loop coj^ = ujk is the spin- 
wave excitation frequency and /"^"^(k, qi, q2) = 
7VA^^(qi,qi -q2)^(k,qi - q2). 



M) 



Here A^^"^ = l/{e^< - 1) ^k = 1 + K'' and A' = N^ 



r(i) 



(1) 



The most significant vertex corrections due to purely 
spin- wave loops are given in Ref. 13 and are those of 
Fig. [191 Their contributions can be calculated at finite 
temperature in a simple way as follows. The contribution 
of the first two diagrams is obtained from the expression 



-9i:ci2 



qi -^" qp 



.'-*' "»'-\ 



gi-ci2 





FIG. 18: The leading vertex corrections (two-loop) to the 
hole Green's function due to phonon (solid wiggly line) and 
spin- wave (dashed- wiggly line) loops. 



given by Eq. [4] by multiplying it with the prefactor (" de- 
fined in Ref. 13 which is the factor that renormalizes the 
spin- wave velocity ^^. The contribution of the third dia- 
gram together with the leading order one loop diagram 
(given by Eq. [Ij) is obtained from the same expression 
given in Eq. [4] by replacing the spin-wave velocity uo^ 
with (1 + C)^/c- The last two- loop diagram of Fig. [19] is 
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k ;^ k-q -; k^ 









( k k':' k-q if k-p >k 



Ju: 



k-q 



FIG. 19: The leading vertex corrections (two-loop) to the hole 
Green's function due to spin- wave loops. 



given by the following expression: 

S(k,u;) = ^p2(q,p)-F(k,q)F(k,p), 
p,q 

F(k,q) = 5(k,q)(iVf G(c. + u;q,k-q) 

+ (l+iV42))G(u;-a;q,k-q)), (11) 



where p2 (k, q) is defined in Ref. fisl. 

We found that the contribution from vertex corrections 
to be small even up to intermediate phonon coupling. 
However, the difference is significant in the strong cou- 
pling limit. Fig. [2Ql shows the spectral function with and 
without the vertex correction using 7 = 0.2t fFig. [2Qlf a)) 
and 0.5t (Fig. EO^b)) and pt = 10. As it can be inferred 
from Fig. [2Qlf b) the lowest energy peak at (^, ^) becomes 
more broadened and there is more rapid transfer of spec- 
tral weight as we move along the (f , f ) — (0, 0) direction. 
Hence, owing to the vertex corrections the "waterfall" - 
like feature, observed in the ARPES experiments, can be 
reproduced using smaller values of the electron-phonon 
coupling. If we include the vertex corrections the (f , f ) 
peak reduces its intensity and at strong enough coupling 
the (0, 0) high energy peak becomes more intense than 
the lowest energy peak. Therefore, we can remain in the 
intermediate coupling regime (e.g., 7 ~ 0.5t shown in 
Fig. [2T]) and still be able to reproduce a broadening simi- 
lar to that observed in ARPES (bottom part of Fig.[2T]). 



VI. CONCLUSIONS 

We have studied the t — J-Holstein model (including its 
t — t' — t" — J extension) at finite temperature within non- 
crossing approximation. We have also included vertex 
corrections due to the electron-phonon coupling. We find 
that the string excitations considered in Ref. 6 to account 
for the water fall- like features of the spectral function ob- 
served in ARPES^i^ are robust even at strong electron- 
phonon coupling and at room temperature. Furthermore, 



1 ' 


1 ' 1 ' 1 ' 1 . 


,/// 


/ 'llO/ /^^^^^^^C'''''''"""--^ "'^''^=^=^^=-. 
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''l\J X^'^^-. ^"^"°"'"'*'---— J 




JJa\ ^ J \ '"*'**---_,_____ 


, , -:jj-^ , , >-_. , 




FIG. 20: Spectral function for (a)7 = 0.2t and (b)7 = 0.5t 
for I3t — 10 and on a 16 x 16 lattice with (dashed lines) and 
without (solid lines) vertex corrections due to the electron- 
phonon interaction. 



the hole spectral function obtained from the NCA treat- 
ment compares well with the reported ARPES intensity 
if we adopt a strong (7/t ~ 1) hole-phonon coupling. 
Namely, it exhibits the same general behavior found in 
Ref. 6, where an artificial spectral broadening was used 
in order to compare with the ARPES data; in the present 
treatment this agreement is achieved without using any 
such artificial broadening procedure. In the calculation 
reported in the present paper, the width and the energy 
dispersion of the lowest energy peak near (7r/2, 7r/2) is re- 
produced and, in addition, we are able to qualitatively re- 
produce the high energy anomaly, i.e., the abrupt down- 
turn in intensity which is characterised by two energy 
scales and the flat featureless intensity between them^. 

Our calculation, where we included the leading-order 
vertex corrections due to the electron-phonon coupling, 
indicates that the vertex corrections are relatively small 
up to an intermediate coupling regime. We also found 
that in the strong coupling limit, they become significant, 
as expected. Furthermore, in order to reproduce the ob- 
served features in the ARPES spectra and intensity plots 
when we included the contribution of the vertex correc- 
tions, the value of the electron-phonon coupling needed 
to achieve the same qualitative agreement was found to 
be smaller than the one needed using our results obtained 
with the non-crossing approximation. This suggests that 
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the qualitative features of the results obtained within the 
non-crossing approximation (which is expected to fail in 
the strong coupling regime) might be valid in the regime 
describing the cuprate materials. 



(71/2,71/2) 



(0,0) 



(jr/2,7i/2) 



FIG. 21: Top: Calculated intensity plot of a 16 x 16 lat- 
tice at P — 10 for 7/t = 0.5 where vertex corrections have 
been included. Bottom: ARPES intensity along the (0, 0) to 
(7r/2,7r/2) direction. 
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